2/3 + 7 from B° — > D t K q -k ± decays at BaBar: a simulation study. 
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We present the results of a simulation study to perform the extraction of 2/3 + 7 from B° — > 
D =F _ft' 7r ± decays through a time-dependent Dalitz analysis of BaBar data. 



I. INTRODUCTION 

A decay of a B° is sensitive to the weak phase 2/3 be- 
cause of the B°B° mixing. If the decay leads to a final 
state which can be reached both through Cabibbo al- 
lowed b — > c (V c b) transitions and Cabibbo suppressed 
b — ► u (V U b) transitions, there is also sensitivity to the 
weak phase 7. Thus we have sensitivity to the total 
phase: 

—V* 2 V 2 V* V 

cd c b cs 

where Vsa are the terms of the CKM quark mixing ma- 
trix [l| 0. 

Up to now the constraints on 2/3 + 7 come from the 
analysis of the two-body decays B° — » D^n and B° — > 
Dp [3|. Another way of measuring 2/3 + 7 is to study 
the B° — > D =F AT°7r ± decays, as proposed in [4j. Fig- 
ure [1] shows the Feynman diagrams for the Cabibbo al- 
lowed and suppressed processes. Only tree-diagram con- 
tributes. The main advantage of these three-body decays 



comes from the possibility of performing a Dalitz anal- 
ysis, which allows to reduce the typical eight-fold ambi- 
guity in the determination of 2/3 + 7 to only a two-fold 
ambiguity Moreover the ratio r of the Cabibbo al- 
lowed over the Cabibbo suppressed amplitudes, which in 
the two-body decays is about 0.02, is expected here to be 
quite large (about 0.4) in some regions of the Dalitz plot. 
This parameter is directly related to the sensitivity to 
2/3 + 7 . The ambition of the analysis technique is to be 
completely model independent, in the sense that at high 
statistics it will not rely on any theoretical assumption: 
all quantities will be directly extracted from the fit on 
data. 

A sensitivity study of the channel can be found in [f| . 
Here, following the same approach, we complete it with 
simulations making use of realistic background distribu- 
tions, resolution and tagging performances taken from 
the BaBar Runl-4 data sample. 



II. DALITZ CONTRIBUTION TO THE TIME 
DEPENDENT EQUATION 




(a) 




The model assumed for the decay parameterizes the 
amplitude A at each point k of the Dalitz plot as a sum 
of two-body decay matrix elements and a non-resonant 
term according to the isobar model Q: 




(2) 

where Ck (ttfc) indicates the Cabibbo allowed (suppressed) 
decay in each point k of the Dalitz plot. Each term of 
the sum is parameterized with an amplitude (oj or a nr ) 
and a phase {5j or 5 nr , where the nr index indicates 
the non resonant). The factor BWl(m,r,s) gives the 
Lorentz invariant expression for the matrix element of a 
(c) (d) resonance j as a function of the position fc in the B Dalitz 

plot; the functional dependence varies with the spin s of 
FIG. 1: Feynman diagrams for the processes contributing the resonance, the mass m and the decay width T. 
to the D*K%± decays; a) B° ^ D-_K*+{K%+) and The ^ dependent evolution can be obtained from 

hiqher K resonances, b) B — > D (D ty )K , c) B — > ,i 1 , • r j.i ci i j- i- 1 j- a 

„ ,,n N _ . 1 r~»o 7^**0/7^+ -,,,0 m, the resolution 01 the behrocdmger equation, leading to 

D s T (_D T K J7r , d) B 1 — » D (D T 7r )K . The processes 




in a,b (c,d) are Vcb (Kit. ) mediated. 



the following expression: 



•Electronic address: polci@lal.in2p3.fr, schunem@lal.in2p3.fr, stocchi@lal.in2p3.fr +eS^ Sm(AiridAt)} , (3) 



2 







r(^) 


a(V cb ) 


a(V ub ) 


S(V cb ) 


8{V ub ) 


D Sj2+ (2573) ± 


2.572 


0.015 


- 


0.02 


- 


70 


D* + (2460) u 


2.461 


0.046 


0.12 


0.048 


30 


30 


D*+(2308) u 


2.308 


0.276 


0.12 


0.048 


70 





/ v | _ ^oyz j 


n sol fifi 


U.UJUO 


± 




n 

u 




lift- (1430)* 
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if* + (143Q) ± 


1.4256 


0.0985 


0.2 









(1680)* 


1.717 


0.322 


0.3 




30 




Non resonant 






0.07 


0.028 





30 



TABLE I: Masses and widths of the resonances considered in 
the model. The last four columns present the chosen values of 
the coefficients aj and 5 in Eq.\^for the Cabibbo allowed and 
Cabibbo suppressed decays respectively. 



with 



C k = 



QK 



A 2 



A 2 



A 2 



2Im(A Ck A Uk e 



i(2/3+ 7 )+r)i(5 c . -S u . ) 



A l k + A u k 



(4) 



where e = +1 (-1) if the tagged initial state is a B° (B°), 
r) = +1 (-1) if the final state contains a D + (D~), and 
At is the proper time interval between the reconstructed 
B (B rec ) and the tagging B (Bt ag ). Tb and Arrid are 
respectively the B° decay width and its mixing frequency. 
Thanks to the presence of the terms BW£ (m, T, s) which 
vary over the Dalitz plot, we can fit the amplitudes (aj) 
and the phases (Sj) of Eq. [21 together with 2/3 + 7 with 
only a two-fold ambiguity. 



III. THE B° -» D^K ^ DALITZ MODEL 

The decay mode under study is assumed to proceed 
through the resonances listed in Table |TJ The K°tt + res- 
onances (K* like) can only come from V c b mediated pro- 
cesses, while both V c b and V u b transitions contribute to 
D~ir + resonances (D** like). Finally the D~K° reso- 
nances (D**-like) come only through V u b mediated pro- 
cesses. Contributions coming from very wide resonances 
like higher K excited states or higher D** excited states 
will be taken into account by a generic "non resonant" 
term. 

The total phase and amplitude are arbitrary, so we 
can chose amplitude unity and phase zero for the mode 
K* + (892) decaying into K a n + . All the other amplitudes 
and phases values are referred to the K* + (892) ones. 
Since the strong phases are not known experimentally, 
their values are chosen arbitrarily, while the values of 
the amplitudes comes from the available measurements 
of related B decay modes and some theoretical consid- 
erations Q. In agreement with the Standard Model we 
assume for the sensitivity study 2/3 + 7 = 2 rad. 

The total number of D t K°tt ± signal events is esti- 
mated using the measured branching fraction and ob- 




m (DKs) 



FIG. 2: Dalitz distribution of the very high statistics Monte- 
Carlo sample of signal events (black dots). Each event is 
weighted with the second derivative with respect to 2f3 + 7 of 
the log-likelihood (colored squares). 



served yield @ 
unit of 100 fb- 



we expect about 250 signal events per 



IV. SENSITIVITY TO 2/3 + 7 

In order to show the regions of the Dalitz plot that 
mostly contribute to the determination of 2/3 + 7, a very 
high statistics Monte-Carlo sample of signal events has 
been generated according to the nominal model described 
in Section HTT1 Since the uncertainty on 2/3 + 7 is: 



1 



02/3+7 — 



J2 8-lnL 



(5) 



9 2 (2/3+7) 



we can weight each event by the quantity weight, the sec- 
ond derivative with respect to 2/3 + 7 °f the log- likelihood 
constructed according to Eq[3j 



weight 



d 2 lnL 
9 2 (2/3 + 7 ) 



(6) 



Note that this likelihood considers only signal contribu- 
tion and does not take into account tagging and resolu- 
tion effects. 

In Figure [21 the weighted distribution expressing the sen- 
sitivity (colored squares) is superimposed to the distri- 
bution of the Monte Carlo events in the m 2 D+K0 versus 
rrvL plane (black dots). The regions with interfer- 
ence between B° -> D**°K° and B° -► D**°K° color 
suppressed processes (diagonal side of the Dalitz) show 
the greater sensitivity to 2/3 + 7. A particularly sensitive 
zone is at the intersection between B° — ► D~K* + and 
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FIG. 3: Average absolute uncertainty on 2/3 + 7 as a function 
of the integrated luminosity. Magenta thick stars refers to the 
model in tabled thin stars are obtained by varying the strong 
phases values; red full circles and green open squares come 
from multiplying the £)**° amplitudes by a factor 1.3 and 0.7 
respectively; red full and green open triangles refer to previous 
models considering r=0.5 and r=0.3 respectively; blue crosses 
refers to the model in table\]\ with no non-resonant component. 
The cyan full squares correspond to the model with ordy non- 
resonant amplitudes for all V u b processes presented in Jj/.- this 
is in contradiction with present data 



the color suppressed B° — ► ]J**°K° (right bottom cor- 
ner of the Dalitz plot). Some sensitivity in the vertical 
left side of the plot is present because of the interfer- 
ence of the Z? S! 2(2573) ± with the excited K resonances 
(^(1430) ± , X 2 *(1430) ± , X*(1680) ± ). 

Since the present measurements (0,@|) are not suffi- 
cient to fix the Dalitz model, the dependence on the de- 
termination of 2(3 + 7 from the assumed Dalitz structure 
of the decay has been evaluated, has shown on Figure [3] 
Here 2/3 + 7 is fit leaving fixed all the other parameters, 
and the impact on the error on 2/3 + 7 of variations in 
the decay model is studied as a function of the luminos- 
ity. The precision on 2/3 + 7 strongly depends on the 
variation of the branching fractions of neutral B into D** 
states within the measured errors. 

An important feature of this method is the possibility 
to have a model independent determination of 2/3 + 7 
by fitting the parameters of the decay model. This is 
shown on Figure [U Considering a sample correspond- 
ing to 1 ab^ 1 of collected statistics, the reference error 
is obtained as previously by fitting only 2/3 + 7 with a U 
the other parameters fixed. Then the other parameters 
which characterize the decay model are progressively re- 
leased. At this point untagged events are useful: even 
if they do not carry information on 2/3 + 7 , still they 
help in determining amplitudes and phases, lowering the 
uncertainty on 2/3 + 7 . 

A first attempt to estimate the effect of the background 



FIG. 4: Average absolute uncertainty on 2/3+7 as a function 
of the fit configuration ("T" stands for tagged and "U" for 
untagged events). The magenta thick stars refers to fit with 
no background while red dots (blue squares) correspond to fits 
where the level of background (fiat in the Dalitz plot) is set at 
50% (30%). 



is done generating it flat over the Dalitz plot. The signal 
over background ratio is fixed to 30% or 50%. The result, 
shown in fig. [H is that the error on 2(3 + 7 increases by 
25% and 50%, respectively. 



V. SIMULATION 

The previous study has been completed performing 
new simulations, which take into account realistic back- 
ground levels and distributions over the Dalitz plot, as 
well as the resolution effects, as obtained from studies on 
BaBar data and Monte Carlo samples. 

Some variables are useful to discriminate between sig- 
nal and background: the beam-energy substituted mass 

niEs = \J (v 7 */^) 2 — Pi? 2 ; the difference between the 
B° candidate's measured energy and the beam energy 
AE = E* B — \fs/2 (the asterisk denotes evaluation in 
the T(45) CM frame); a fisher discriminant FI combi- 
nating some event shape variables; the mass mo of the 
D meson. We define as the product of the probabil- 
ity density functions (PDFs) of these variables for each 
event i and for each component j: signal (Sig), contin- 
uum background (Cont), combinatoric BB decays (BB 
background) and BB events that peak in uies but not 
in AE signal region (denoted peaking BB background: 
Peak). 

The complete expression of the likelihood function, 
taking into account also the time dependence and the 
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a(V ub ) 


S(V cb ) 


S(V ub ) 


D s , 2(2573)* 
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0.002 ± 0.011 
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± 155 


D5(2460) u 


0.11 ± 0.01 


0.046 ± 0.019 
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DS(2308) U 


0.13 ± 0.02 


0.047 ± 0.023 


70 ± 11 


38 ± 31 


K* (H<n\^ 


i 
± 




n 
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K£ ( 1430) * 


0.60 ± 0.02 




81 ± 2 




#5(1430)* 


0.20 ± 0.03 




± 2 




#■•(1680)* 


0.30 ± 0.01 




28 ± 4 




Non resonant 


0.09 ± 0.02 


0.070 ± 0.030 


357 ± 13 


48 ± 24 



TABLE II: Results of the fit simulation corresponding to an 
integrated luminosity of 10ab~ x 



Dalitz distribution is: 



B° tag 



B° tag 



(7) 



where N c — 7 is the number of tagging categories, 

is the likelihood function for an event in the tag- 
ging category i with B tag = B° (B tag = B°) and: 



K § V BS Tl BB y BB (8) 

+N % Cont Vcont'r±,ContyCont (9) 
+ N Peak D PeakT± i PeakyPeak 

Here V is the time-dependent Dalitz PDF for signal 
(eq. [3]), T> and T indicate the Dalitz and the time pa- 
rameterizations for the backgrounds. 

The simulations are performed taking the T>, T and y 
shapes for the backgrounds from BaBar data and Monte 
Carlo simulations, while for the signal V the model previ- 
ously described is taken into account. About 1000 sam- 
ples of signal and background are generated and fitted 
according to these PDFs. 

The results show that an integrated luminosity of 
354 (corresponding to Runl-4 BaBar data) is not 



enough to extract all amplitudes and phases from the fit. 
In particular the non resonant component is not correctly 
determined. 

For this reason other simulations are performed intro- 
ducing a parameter r = 0.3, so that for each resonance 
the relation: A(V u h) = r x A(V c b) holds. The choice 
of the 0.3 value for r is suggested by the result of the 
B° -> D°K*° analysis [|, giving a limit r < 0.4 @ 90% 
CL. 

The simulation at 354 fb' 1 in this configuration works 
if the amplitudes and phases for the D** and the non 
resonant component are fixed to the model value, as well 
as all the V u h phases. In this case the V c b amplitudes 
are determined with about 30% error, and the error on 
2/3 + 7 is 01 about 75 degrees. Some biases are observed 
in the amplitudes determination, which are consistently 
reduced if the simulation is performed at a statistics of 
lab^ 1 : in this case the error on 2/3 + 7 is about 57 
degrees. 

The last test performed is a simulation at very high 
statistic (10a6 _1 ), where we do not use the parameter 
r, but we fit directly all amplitudes and phases gener- 
ated at the values presented in table |TJ Results, shown 
in table HT1 show that all the parameters can be fitted, so 
that 2/3 + 7 is determined with an error of about 14 de- 
grees in a completely model independent way. Of course 
the estimation on the error on 2/3 + 7 depends strongly 
on the Dalitz model assumed: reality could be different, 
depending on the relative abundances of the resonances. 

In conclusion, a model independent analysis method 
for the determination of 2/3 + 7 has been presented 
and validated using realistic simulation, based on the ob- 
served shapes of the discriminating variables on BaBar 
data and simulated samples as well as on BaBar tagging 
and resolution performances. At the luminosity of la6 _1 
,2/3 + 7 can be determined with an error of about 57 
degrees. 
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